Neutron skin thickness in droplet model with surface width dependence: 
indications of softness of the nuclear symmetry energy 
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We analyze the neutron skin thickness in finite nuclei with the droplet model and effective nuclear 
interactions. The ratio of the bulk symmetry energy J to the so-called surface stiffness coefficient 
Q has in the droplet model a prominent role in driving the size of neutron skins. We present a 
correlation between the density derivative of the nuclear symmetry energy at saturation and the 
J/Q ratio. We emphasize the role of the surface widths of the neutron and proton density profiles in 
the calculation of the neutron skin thickness when one uses realistic mean-field effective interactions. 
Next, taking as experimental baseline the neutron skin sizes measured in 26 antiprotonic atoms along 
the mass table, we explore constraints arising from neutron skins on the value of the J/Q ratio. 
The results favor a relatively soft symmetry energy at subsaturation densities. Our predictions are 
compared with the recent constraints derived from other experimental observables. Though the 
various extractions predict different ranges of values, one finds a narrow window L ~ 45-75 MeV for 
the coefficient L that characterizes the density derivative of the symmetry energy which is compatible 
with all the different empirical indications. 
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I. INTRODUCTION 

Neutron skin thickness is the name in common usage 
to refer to the difference between the root-mean-square 
(rms) radii of the neutron and proton density distribu- 
tions of atomic nuclei: 

Ai?„, = (r^)y^ - {rYj\ (1) 

Experimentally, the value of the proton rms radius 
(r^)^^ is obtained from the charge radius. The latter 
has been measured by electron-nucleus elastic scattering 
with high accuracy (often the accuracy in charge radii 
is better than 1% [1]). In contrast, our knowledge of 
the neutron distribution in nuclei and of its rms radius 
{r'^)r/'^, as well as our knowledge of Ai?„p, is till date less 
precise. This situation looks inadequate in anticipation 
of the next generation of rare ion accelerator facilities 
that are planned for the synthesis and study of exotic 
nuclei, as in RIKEN centers (Japan), in FAIR at GSI 
(Germany), in the GSR at HIRFL (Ghina), or in FRIB 
at MSU (USA) . Without securing sufficient knowledge of 
the neutron distribution in stable nuclei, the prospects of 
nuclear structure theory in this thriving domain may be 
compromised. It is expected that parity-violating elec- 
tron scattering will provide in the nearby future a leap 
forward in the quest for high precision determinations of 
the neutron radius in heavy nuclei [2] . 
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The calibration of the neutron skin thickness of nu- 
clei also is one of the problems at the forefront of nu- 
clear structure by reason of being intimately correlated 
with the nuclear symmetry energy. Indeed, the symme- 
try energy is a fundamental quantity in nuclear physics 
and astrophysics, because it governs at the same time 
important properties of very small entities like atomic 
nuclei and of very large objects like neutron stars [3]. 
One of the crucial properties of the symmetry energy, 
which still is not sufficiently well constrained, is its de- 
pendence on the nuclear density. It is relevant in the 
astrophysical context in order to understand a wealth of 
phenomena [3~5], including supernova explosions, neu- 
trino emission, and the cooling mechanism of protoneu- 
tron stars, as well as mass-radius relations in neutron 
stars. Moreover, the density content of the symmetry en- 
ergy eventually relates to basic issues in physics. This is 
the case of precision tests of the Standard Model through 
atomic parity non-conservation observables [6] , and even 
of studies on constraining a possible time variation of 
the gravitational constant [7]. In terrestrial laboratories, 
the available tools to delineate the density dependence 
of the symmetry energy at saturation and subsatura- 
tion densities, include the interaction potential between 
neutron-rich nuclei [8], observables like isospin diffusion 
and isoscaling in heavy-ion reactions at intermediate en- 
ergies [9-21], different modes of collective excitations of 
nuclei [22-26], and, of course, data on the binding and 
structure of neutron-rich nuclei and on their neutron skin 
thickness. 

In the literature there exist several theoretical formula- 
tions to investigate the neutron skin thickness of neutron- 
rich nuclei and its connections with the symmetry energy. 
This is the case, for instance, of methods based on the 



2 



droplet model [27, 28], on the concept of surface symme- 
try energy [5, 29, 30], thermodynamical arguments [31], 
rmclconic density form factors [32], mean- field analyses 
[33-35] , or studies in the spirit of the Landau-Migdal ap- 
proximation [36]. It has been shown that the neutron 
skin thickness in heavy nuclei, like ^"^Pb, calculated in 
mean-field models with either non-relativistic or relativis- 
tic effective nuclear interactions, displays a linear corre- 
lation with the slope of the neutron equation of state 
(EOS) obtained with the same interactions at a neutron 
density p « 0.10 £m~^ [37-39]. A similar correlation ex- 
ists between Ai?„p and the density derivative of the bulk 
symmetry energy [13-16, 33, 40, 41], as the latter is a 
measure of the pressure difference between neutrons and 
protons. These correlations have been exploited in recent 
years to gain a deeper understanding of the isospin prop- 
erties of the effective nuclear interaction and to relate 
them with nuclear and astrophysical observations. 

The rms radius of neutron densities in nuclei has been 
measured with hadronic probes such as proton-nucleus 
elastic scattering [42-45] or inelastic scattering excitation 
of the giant dipole and spin-dipole resonances [46, 47]. 
On the other hand, antiprotonic atoms are helpful to 
probe the size of the neutron skin of nuclei from the 
fact that the nuclear periphery is very sensitive to an- 
tiprotons in the normally electronic shell. Experimen- 
talists combine two different techniques in this case [48- 
50], namely, the measurement of the antiprotonic X-rays 
which determine the atomic level shifts and widths due 
to the strong interaction, and the radiochemical analy- 
sis of the yields after the antiproton annihilation. The 
values of the neutron skin thickness of 26 stable nuclei 
from ■^^Ca to ^38^ deduced from antiprotonic atoms data 
by Trzcinska et al. [48, 49] follow a roughly linear trend 
with the overall relative neutron excess I = [N — Z)/A 
of these nuclei. This trend can be fitted by the rela- 
tionship Ai?„p = (0.90 ± 0.15)/ + (-0.03 ± 0.02) fm as 
discussed in Refs. [48, 49]. As mentioned, all neutron skin 
thickness measurements have relatively large uncertain- 
ties in comparison with charge radii, and sometimes the 
results from different experimental techniques are not to- 
tally consistent among them [47, 49]. The neutron skin 
sizes determined in Refs. [48, 49] from the analysis of 
antiprotonic atoms are till date the largest set of uni- 
formly measured values of Ai?„p all over the periodic 
table (40 < A < 238). Due to this reason, we shall use 
hereinafter these data as the experimental benchmark for 
our calculations. 

The droplet model (DM) describes in a physically 
transparent way nuclear radii and relates them directly 
with basic properties of the nuclear interactions. In the 
present paper we study the neutron skin thickness of 
atomic nuclei with the DM using various effective nu- 
clear interactions of the Skyrme, Gogny, and relativistic 
mean-field (RMF) type. The present work extends with 
a new analysis and perspective a first presentation of our 
study made in Ref. [51]. Here, we will show that the ratio 
of the DM parameters J and Q, which drives the value 



of the neutron skin thickness in heavy nuclei, is corre- 
lated with the slopes in density of the nuclear symmetry 
energy and of the EOS of neutron matter. We compare 
the DM values for the neutron skin thickness with the re- 
sults obtained in self-consistent extended Thomas-Fermi 
(ETF) calculations of finite nuclei [54, 58, 62, 63], since 
both methods are free of shell effects. A non-negligible 
role of the contribution of the difference in the surface 
widths of the neutron and proton density profiles is no- 
ticed. Next, we use the experimental neutron skin thick- 
ness measured in antiprotonic atoms to explore the range 
of possible values of the ratio J/Q that are favored by 
neutron skins. With these values we can predict some 
properties of the density dependence of the nuclear sym- 
metry energy. Our results are compared with the con- 
straints recently obtained in the literature using other 
observables and methods. 

The present paper is arranged as follows. In the sec- 
ond section we study the neutron skin thickness of heavy 
nuclei on the basis of the DM [27, 52, 53] and show a cor- 
relation that links the value of the ratio J/Q, which gov- 
erns the neutron skin thickness of nuclei, with the slope 
of the symmetry energy in bulk matter at saturation. In 
the third section, the contribution of the surface widths of 
the neutron and proton density distributions to the neu- 
tron skin thickness is analyzed with the DM using non- 
relativistic and covariant mean-field nuclear interactions. 
In the fourth section we estimate possible constraints on 
the density dependence of the nuclear symmetry energy 
on the basis of the DM and the experimental data on 
neutron skin sizes derived from antiprotonic atoms. We 
discuss the present results in comparison with the recent 
constraints obtained from various observables and meth- 
ods. Finally, the summary and our conclusions arc laid 
in the fifth section. We outhne the procedure for the 
calculation of the Q coefficient in the Appendix. 



II. THE FRAMEWORK 

A. Neutron skin thickness in the droplet model 

In the DM of average nuclear properties [52-54] the 
neutron skin thickness of a finite nucleus is computed 
from the expression [27, 52] 



■np 



(2) 



where e^Z/lQJ is a correction due to the Coulomb in- 
teraction, R = tqA^^^ is the nuclear radius, and 6„ and 
bp are the surface widths of the neutron and proton den- 
sity profiles. In the "standard" version of the DM it is 
assumed that bn = bp — I fm [27, 28, 52], which implies 
a vanishing surface width correction to the neutron skin 
thickness. 

The quantity t in Eq. (2) represents the distance be- 
tween the neutron and proton mean surface locations. 
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FIG. 1: (Color online) For several nuclear mean-field models, 
existing correlation between the neutron skin thickness Ai?„p 
in ■^'^^Pb and the ratio J/Q (left panel), and between AR„p in 
^°*Pb and the slope of the symmetry energy L (middle panel). 
The correlation between the coefficient L and the ratio J/Q 
is also shown (right panel). In the present figure, ARnp has 
been computed with Eq. (1) from the rms radii of quantal self- 
consistent calculations for the indicated mean-field models. 



This distance is computed as [27, 52] 

r 1/3 

3 J vn 
* = 9*^0 n TT 

4(3 
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where I — {N — Z)/A, J is the symmetry energy co- 
efficient at saturation, Q is the surface stiffness coeffi- 
cient, and ci — 3e^/5ro. The coefficient J represents 
with a very good accuracy the energy cost per nucleon 
to convert all protons into neutrons in symmetric infi- 
nite nuclear matter at saturation density po- The sur- 
face stiffness coefficient Q measures the resistance of the 
system against separation of neutrons from protons to 
form a neutron skin. To extract Q from an effective nu- 
clear interaction requires to perform calculations of asym- 
metric semi- infinite nuclear matter (ASINM). Therefore, 
the calculated value of Q may depend somewhat on the 
type of approach, such as the Hartree-Fock or Thomas- 
Fermi methods, employed to describe the nuclear surface 
[30, 54-58]. 

From Eq. (3), one sees that the leading contribution 
to t in large nuclei is the term ^ro{J/Q)I. Thus, the 
DM suggests that one can expect a correlation between 
ARnp and J/Q in heavy nuclei. We illustrate this fact in 
the left panel of Fig. 1. We depict there the neutron skin 
thickness of ^"^Pb obtained from self-consistent quantal 
calculations with the Skyrme and Gogny Hartree-Fock 
methods as well as with the RMF Hartree approach. The 
results are shown as a function of the value of the J/Q 
ratio for various mean-field effective interactions. The 
J values of the selected interactions (about 27-32 MeV 
in the non-relativistic forces and about 32-45 MeV in 
the covariant forces, see Table I) cover widely the plau- 
sible physical range of the bulk symmetry energy. The 
values of Q used in this work have been extracted from 



TABLE I: Saturation density po and J, L, Ksym and Q pa- 
rameters of the Skyrme and Gogny forces as well as RMF 
parameter sets used in this work. 
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ASINM calculations performed in the extended Thomas- 
Fermi (ETF) approach as described in the Appendix (see 
also Ref. [58]). Even if the shell effects, present in the 
mean- field calculations of ARnp, are not built in the DM 
[27], one observes a considerably linear correlation be- 
tween the values of ARnp and J/Q. It should be pointed 
out that while all the effective interactions have been 
accurately calibrated to data on binding energies and 
charge radii, and describe these properties very success- 
fully, they predict widely different values of ARnp, as we 
see in the present case of ^"^Pb. This underlines the fact 
that the isospin sector of the effective interactions is little 
constrained. 



B. Properties of the nuclear symmetry energy 

Let us consider the energy per particle e(p, 6) in asym- 
metric infinite nuclear matter of total density p ^ Pn + Pp 
and relative neutron excess 6 — {pn — Pp)/p, where pn 
and Pp stand for the neutron and proton densities, re- 
spectively. The general expression 



e(p, 6) = e(p, ,5 = 0) + Csymip)S^ + 0(6^) 



(4) 



defines the symmetry energy coefficient Csym (p) of a nu- 
clear EOS at the density p. This expression is particu- 
larly useful because Csym (p) dominates the corrections to 
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the symmetric limit for all values of S, especially at the 
subsaturation densities of relevance for finite nuclei [59]. 
Actually, Cgymip) provides with excellent accuracy the 
difference between the binding energies of pure neutron 
matter {S = 1) and symmetric matter (5 = 0). 

It is customary, and insightful, to characterize the hc- 
havior of an EOS around the saturation density po by 
means of a few bulk parameters calculated at the satu- 
ration point, as in the formula [10, 13-16, 52, 59-61] 



and of the symmetry energy 



i{p,d) « + — e + 



J-Le 



K. 



sym ^2 



(5) 



where e = {po — p)/^Po expresses the relative density 

displacement from pQ. Here, the quantities a^, and Ky 
denote the energy per particle and the incompressibil- 
ity modulus of symmetric nuclear matter. One has 



i(Po) 



J. The DM coefficients L and Ks 



respectively, proportional to the slope and the curvature 
of the symmetry energy coefficient Csym{p) at saturation 
density: 



3po- 
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sym 
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ip) 
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The quadratic expansion Csym 

in Eq. (5) is often a reliable representation of the actual 
value of the Cgymip) coefficient at densities roughly be- 
tween po/2 and 2po [59]. For instance, in the case of a 
typical subsaturation density value p = 0.10 fm~^, one 
finds that the above quadratic expansion of Csym{p) dif- 
fers from the exact Csym{p) by less than 1% in many dif- 
ferent nuclear mean field forces [51]. These facts point to 
the usefulness of investigating parameters such as L and 
Ksym for the characterization of the density dependence 
of the symmetry energy. 

The values of the DM coefficients J, L, and Ksym 
for the non-relativistic forces and the RMF parameter 
sets considered in this work are given in Table I. In the 
Skyrme and Gogny effective interactions the symmetry 
energy coefficient at saturation J takes values around 30 
MeV. The RMF parameterizations have larger values of 
J, also with a larger spread. The slope (L) and the curva- 
ture {Ksym) of the symmetry energy at saturation take 
even more widely scattered values among the different 
interactions. The consequence is that all mentioned nu- 
clear models predict a different behavior of the symmetry 
energy at subsaturation densities, what can be seen e.g. 
in Fig. 1 of Ref. [13]. 

As is known, the density dependence of the symme- 
try energy near saturation tends to be much softer in 
the non-relativistic forces than in the covariant meson- 
exchange models of nuclear structure (see the values of 
L in Table I). In particular, taking into account Eq. (5), 
the slope of the neutron EOS 



dCsym^P) 

dp 



L 



K. 



sym 



3po 
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calculated at densities p close to the saturation value po 
differ considerably between models. We also know that 
Ai?„p in 2'^^Pb shows a linear dependence with these 
slopes at some subsaturation density p ~ 0.10 fm~"^ [37- 
40]. Therefore it is reasonable that the neutron skin 
thickness in ^''^Pb and the leading term L/ipo of Eqs. 

(7) and (8) are related. As far as po does not change 
much in the different effective forces, a correlation be- 
tween ARnp in ^''^Pb and the DM coefficient L is thus 
expected [13-15, 33, 41] and we display it in the middle 
panel of Fig. 1. 

From the discussed results of ARnp versus J/Q and 
of ARnp versus L in Fig. 1, a correlation between the 
DM coefficient L and the J/Q ratio is to be expected 
too. Note that L rules the density dependence of the 
symmetry energy of the nuclear equation of state [Eq. 

(8) ], and that Q governs the thickness of the neutron skin 
of finite nuclei [Eqs. (2) and (3)]. The correlation between 
L and J/Q can be seen in the right panel of Fig. 1. This 
correlation shows that the value of the slope L of the 
symmetry energy at the saturation density increases with 
the value of the ratio between the bulk symmetry energy 
coefficient J and the surface stiffness coefficient Q. The 
trend is considerably linear among the various nuclear 
effective interactions. 

Previous literature [56, 58] has shown that the system- 
atics of experimental binding energies relates increasing 
values of J with decreasing values of Q in nuclear effec- 
tive interactions whose parameters have been adjusted to 
describe experimental data. We note this same trend in 
Table I, where the RMF sets that in general have larger J 
values also tend to have smaller Q values than their non- 
relativistic counterparts. A smaller Q coefficient means 
that it is easier to develop a neutron skin in finite nu- 
clei. Consistently, the neutron skin thickness in ^^^Pb 
(or any other heavy nucleus) is usually larger when com- 
puted with a RMF parameter set than when computed 
with a non-relativistic force. These facts, and the noticed 
correlation between L and J/Q, allow one to interpret in 
a qualitative way within the DM, the correlation pointed 
out in Refs. [37-40] between the slope of the symmetry 
energy at some subsaturation density p ~ 0.10 fm~^ and 
the neutron skin thickness in a heavy nucleus. In a re- 
cent work [51] we have investigated further the relations 
of the neutron skin thickness with the parameters L and 
Ksym that characterize the density dependence of the 
symmetry energy around saturation. 



III. SURFACE WIDTH CONTRIBUTION TO 
THE NEUTRON SKIN THICKNESS 



de{p, 6 = 1) 
dp 



L 

3po 



Ky + Ks 
3po 



(7) 



The neutron skin thickness values derived from mea- 
surements performed in antiprotonic atoms have been ob- 



5 



tained in Refs. [48, 49]. It is assumed that the neutron 
skin is due to an enhancement of the neutron surface 
width with respect to the proton surface width, and that 
the mean location of the proton and neutron surfaces in 
these nuclei are the same. This situation corresponds 
to the so-called "neutron halo- type" distribution [48] . It 
has been shown that the same set of experimental values 
of neutron skin thickness can be explained with similar 
quality as in [48] by means of the "standard" version of 
the DM (where 5„ ~ bp) [28]. The latter case assumes 
that the peripheral neutrons are concentrated at the neu- 
tron surface, which is shifted with respect to the proton 
surface, and that both the neutron and proton density 
distributions have the same surface width. This is rather 
the pattern of the so-called "neutron skin-type" distribu- 
tion according to Ref . [48] . 

The analysis of neutron and proton densities calcu- 
lated with nuclear mean field interactions carried out in 
Ref. [32] by means of the Helm model, points out that 
self-consistent mean field densities show a mixed charac- 
ter between the "neutron halo" and "neutron skin" pat- 
terns. This means that, actually, the self-consistent neu- 
tron and proton density profiles obtained with nuclear 
effective interactions differ not only in the mean location 
of their surfaces but also in their surface widths. In the 
following we shall see that similar conclusions are found 
from the calculations of the neutron skin thickness per- 
formed in the DM with formula (2). It will turn out that 
the surface width contribution oc {bf^ — b^) in the DM ex- 
pression (2) for the neutron skin thickness, which arises 
from bn 7^ 5p, is necessary to reproduce the neutron skin 
thickness values calculated from the definition (1) using 
self-consistent densities of finite nuclei obtained with the 
ETF approach in mean-field theory, for non-relativistic 
forces as well as for relativistic parameterizations. 

In Fig. 2 we display by empty symbols, as a function 
of the overall relative neutron excess / = {N — Z)/A, 
the neutron skin thickness predicted by the "standard" 
version of the DM (namely, Eq. (2) with 6„ = bp) using 
some well-known effective forces. The nuclei are those 
from ^°Ca to ^38^ measured in the experiments with an- 
tiprotonic atoms [48, 49], and which were studied with 
the DM in Ref. [28]. The values shown in Fig. 2 have 
been computed using the SIII and SkM* Skyrme forces 
and the NL-SH and NL3 RMF parameter sets, as suit- 
able examples. We have chosen these four parameter sets 
for display in Fig. 2 because they span the whole range 
of values of the ratio J/Q of nuclear interactions that 
describe reasonably well the ground-state properties of 
finite nuclei, having a bulk symmetry energy coefficient 
J between 28 MeV and 37 MeV (see Fig. 1 and Table I). 
The DM results for Ai?„p obtained with the other mean 
field interactions considered in this work that have a J 
coefficient between the values of SIII and NL3, also lie 
within the window of results delimited by the SIII and 
NL3 interactions in Fig. 2. 

The values of Ai?„j, predicted by the DM are com- 
pared in Fig. 2 with the values that we obtain from 
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FIG. 2; (Color online) The neutron skin thickness predicted 
by the "standard" version of the DM (Eq. (2) with bn = 
bp) is compared with the result obtained from self-consistent 
ETF calculations of finite nuclei, in four illustrative mean-field 
parameter sets. The nuclei considered are those investigated 
experimentally with antiprotonic atoms in Refs. [48, 49] and 
have masses 40 < yl < 238. 



self-consistent ETF calculations in finite nuclei (filled 
symbols). Both models do not incorporate shell effects. 
In Fig. 2 we have used the ETF approach in the non- 
relativistic [62] and relativistic [63] frameworks to com- 
pute ARnp- We have calculated these ETF values of 
ARnp by application of Eq. (1) with the rms radii of the 
self-consistent neutron and proton densities obtained in 
each isotope. From Fig. 2 two significant points stem. 
First, the predictions of the DM in the "standard" form 
{bn ~ bp) systematically undershoot the ETF neutron 
skin thickness computed in finite nuclei with the selected 
effective nuclear interactions. In particular, this trend 
is reinforced with growing neutron excess /. Second, it 
can be observed that for a given nucleus the difference 
between the ETF value of ARnp computed with (1) and 
the value provided by the "standard" DM prescription 
is slightly larger in the RMF parameter sets than in the 
Skyrme forces. Altogether, these facts suggest that in the 
mean-field interactions the surface width contribution to 
the DM formula for ARnp does not vanish, and that this 
contribution has some dependence on the ratio J/Q of 
the force. 

To apply the full Eq. (2) to compute the neutron 
skin thickness including the surface width correction, 
one needs to evaluate the neutron and proton surface 
widths in finite nuclei. In practice, there is not conclu- 
sive experimental evidence on the difference of the sur- 
face widths bn and bp of the nuclear density distributions 
[27, 28, 32, 50, 54]. To estimate bl - 5^ we will there- 
fore rely on theoretical guidance as a surrogate. Within 
the context related to the DM and the leptodermous ex- 
pansion of a finite nucleus [52-54], the surface properties 
in finite nuclei can be extracted from ASINM calcula- 
tions. The semi-infinite geometry does not include shell. 
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Coulomb or finite-size efi^ects. We will use here the ETF 
method including fi^ corrections for describing ASINM, 
since the ETF method is free of the Friedel oscillations of 
the quantal densities [30, 57]. We summarize in the Ap- 
pendix the basic aspects of the procedure. More details 
can be found in Ref. [58]. From Eqs. (14)-(17) of the Ap- 
pendix, we can obtain the values of the surface widths 6„ 
and bp in ASINM with a given relative neutron excess in 
the bulk Sq. In the DM, these surface widths correspond 
to the values 5„ and bp in finite nuclei if Sq is calculated 
from the overall relative neutron excess / of the nucleus 
through the following relation [52, 54, 58]: 



1 



9 J 



^-1/3 



(9) 



which takes into account the Coulomb correction. 

Once the neutron and proton surface widths in finite 
nuclei are known, we can compute their contribution to 
the neutron skin thickness, which reads [cf. Eq. (2)] 



3_5_ 
5 2^ 



(10) 
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FIG. 3: (Color online) Upper panel: the same as in Fig. (2) 
but here the DM values include a non- vanishing surface width 
contribution AR!^ [Eq. (10)] with b„ and bp obtained from 
ASINM calculations as described in the text. Lower panel: 
the surface width contribution Ai?^p (the vertical scale pro- 
portionality is the same as in the upper panel). 



The corresponding values of Ai?f™ for the nuclei con- 
sidered in Fig. 2 are displayed in the bottom panel of 
Fig. 3. It is worth to point out that the calculated Ai?™ 
values show, for each nuclear interaction, a well-defined 
increasing linear trend as a function of the overall relative 
neutron excess / of the nuclei. 

The neutron skin thickness predictions of the DM when 
one includes the surface width contribution (2) are dis- 
played in the top panel of Fig. 3 by empty symbols. Note 
that the results correspond to adding the values AR^ 
shown in the bottom panel of this figure to the DM values 
that we have displayed in Fig. 2. As done in Fig. 2, we 
compare in the top panel of Fig. 3 the DM results with 
the self-consistent ETF calculations of Ai?„p [62, 63]. 
One now observes an improved and remarkable agree- 
ment between the DM predictions and the self-consistent 
ETF values computed with the same interaction, stem- 
ming from the inclusion of the calculated AiJ^^ contri- 
bution. It is interesting to note that the neutron skin 
thickness obtained with Eq. (1) from the rms radii of the 
self-consistent ETF calculations in finite nuclei shows a 
well defined increasing linear tendency with the relative 
neutron excess /, similarly to the case of the results of 
the DM and in consonance with the trend of the experi- 
mental values derived from antiprotonic atoms [48]. 

The lower panel of Fig. 3 also suggests that, for a given 
nucleus, Ai?f™ grows when the J/Q ratio of the nuclear 
interaction increases (see Table I) . To analyze this behav- 
ior in more detail, we fit AR^ by means of a law a^^I, 
which defines the slope cr""' of Ai?™ with respect to the 
relative neutron excess /. This slope is displayed in Fig. 4 
as a function of the J/Q ratio for different interactions. 
The slopes cr*™ lie inside a band limited by two straight 




J/Q 



FIG. 4: (Color online) Average slope of Ai?^p with respect 
to / for various nuclear mean-field models as a function of 
J/Q. All the data lie in the area limited by the marked lines 
(7="' = 0.3J/Q + 0.07 fm and cr="' = 0.3J/O - 0.05. 



lines, corresponding to the equations (t"™ — 0.3J/Q+0.07 
fm (left) and a'"" = 0.3J/Q-0.05 fm (right). It is worth 
to notice that all considered Skyrme forces have slopes 
cr*"' below 0.25 fm, whereas the analyzed RMF models 
have slopes cr'"" always above this value. 

One relevant conclusion is that bulk and finite nuclei 
properties described through successful theoretical mean- 
field models constrain the possible values of the surface 
width contribution to the neutron skin thickness between 
the limits portrayed in Fig. 4. From this figure it can 
be deduced that the surface width contribution to the 
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neutron skin thickness AR,f^ has, on top of a global in- 
creasing trend with J/Q, a more involved dependence on 
the parameters of the effective nuclear interactions. For 
instance, on the one hand, the RMF force FSUGold and 
the Skyrme force SkM* have almost the same J/Q ratio 
(see Table I). However, the predicted values of a*"' are 
clearly different for both interactions as can be seen in 
Fig 4. On the other hand, it is possible to find different 
interactions which have almost the same slope cr"™, and 
therefore the same Ai?^™, but with different values of the 
J/Q ratio. Some examples of this fact in Fig. 4 are the 
Skyrme forces SIII and SGII (with slopes of 0.16 fm and 
0.15 fm, respectively), and the RMF parametcrizations 
FSUGold and NL3 (with slopes 0.30 fm and 0.31 fm). 
These facts suggest that it is possible to find the same 
total neutron skin thickness by combining a small value 
of the J/ Q ratio in t [see Eq. (3)] with a large Ai?*™ con- 
tribution or, vice versa, by combining a large J/Q ratio 
in t with a small Ai?^^ contribution. 



IV. ESTIMATES OF THE DENSITY CONTENT 
OF THE SYMMETRY ENERGY 

As we have pointed out, the nuclear symmetry energy 
and the properties of the EOS of neutron-ric;li matter 
are of increasing importance in both nuclear physics and 
astrophysics. There is a significant recent effort in the 
community towards constraining the values of the pa- 
rameters that characterize the density dependence of the 
symmetry energy in the subsaturation regime of the EOS. 
The new developments come in particiilar from the inves- 
tigation of isospin-sensitive observables in intermediate- 
energy heavy-ion collisions [9-21] and in nuclear reso- 
nances [22-25]. Obviously, the different studies do not 
deal with exactly the same regimes of density and energy. 
One also has to keep in mind that the connection of the 
experiments with the EOS is not at all trivial; it often 
requires of extrapolations of the measured data, which 
imply a model dependence. Therefore, it is important to 
further investigate indications from those and other ex- 
perimental probes of the symmetry energy, with different 
methodologies, as well as to study the interplay between 
the constraints derived from the different analyses. 

In the present section we want to apply the experience 
gained in the DM study of the neutron skin thickness 
performed in the previous sections, to estimate possible 
constraints on the density dependence of the nuclear sym- 
metry energy as suggested by neutron skin data. To this 
end, we shall first obtain the range of values of the J/ Q 
ratio which are compatible with the neutron skin thick- 
ness derived from the experimental data in antiprotonic 
atoms. We recall that this set of data for 26 stable nuclei 
is till date the largest set of uniformly measured neutron 
skins spanning the mass table. Once the estimated range 
of values for the ratio J/ Q will be found, the constraints 
derived from the neutron skin data on the DM parame- 
ter L, related to the density derivative of the symmetry 



energy, will be determined from the existing linear corre- 
lation between the values of L and J/Q. This correlation 
has been displayed in Fig. 1 and is a general feature of 
mean field interactions that have been adjusted to re- 
produce with good accuracy binding energies and charge 
radii (often among other properties) of nuclei across the 
periodic table. 

A. Constraints on the J/Q ratio 

From the previous section we know that the surface 
width part in Eq. (2) gives a non-negligible contribu- 
tion to neutron skins in effective nuclear interactions. 
This contribution is needed to reproduce the neutron skin 
thickness values computed self-consistently in ETF calcu- 
lations of finite nuclei. We have also seen that to leading 
order, both the mean location of the neutron and pro- 
ton surfaces (3) and the surface width correction (10) 
are basically driven by the value of the J/Q ratio. The 
discussions in the previous sections suggest to fit the ex- 
perimental AJ?^pP data by means of the following DM 
inspired ansatz: 

where t is given by Eq. (3). The second term is the surface 
width contribution. It is parameterized to reproduce the 
dashed lines on Fig. 4, with c = 0.07 fm or c = —0.05 fm. 

With the ansatz (11), we will use J/Q as an open pa- 
rameter. It will be constrained by a least-squares min- 
imization from the experimental values ARj/jJ' derived 
from the analysis of antiprotonic atoms [48, 49]. We 
note that Eq. (11), as well as t given by Eq. (3), de- 
pends on the particular values of the symmetry energy 
at saturation J and of the nuclear matter radius tq. We 
fix these quantities to the empirical values J = 31.6 MeV 
and tq = 1.143 fm (the latter corresponds to a satura- 
tion density po = 0.16 fm""^). We consider the values 
c = 0.07 fm and c — —0.05 fm in Eq. (11), discussed in 
connection with Fig. 4, to simulate the upper and lower 
bounds of the window of the theoretical predictions for 
cr*"' obtained according to mean-field models of nuclear 
structure. In the x^-minimization we have weighted each 
Ai?JjpP — ARnp difference by the inverse of the associated 
experimental uncertainties. That is, in practice we have 
minimized the quantity 

/A^^i)-A^^y 



where ARnp is calculated with Eq. (11) and the £,i denote 
the imcertainties of the experimental data. 

The fits to experiment give J/Q = 0.667 ± 0.047 with 
c = 0.07 fm and J/Q = 0.791 ± 0.049 with c = -0.05 
fm (i.e., a range 0.62 < J/Q < 0.84). The quoted 
uncertainties in the J/Q predictions correspond to the 
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FIG. 5: (Color online) The values of ARnp obtained from 
Eq. (11) by fitting the J/Q ratio, using c — 0.07 fm (circles) 
and c — —0.05 fm (squares), to reproduce the Ai?^p'' values 
measured in antiprotonic atoms (dots with error-bars). The 
average value of AJJ^p*' is marked by the dotted line. 



value of one standard deviation associated with the fit 
made through Eq. (12). To check our method of mini- 
mization and error estimation, we have apphed the same 
procedure to make a linear fit ml -I- n of the exper- 
imental data AW^pP. In this case we have obtained 
ARnp = (0.901 ± 0.147)/ -I- (-0.034 ± 0.023) fm, which 
fully agrees with the result quoted by the experimental- 
ists [48, 49]. 

The results for ARnp from our fits compared with the 
experimental data AR'fj^P are displayed as a function of 
/ in Fig. 5. Both extractions of J/Q, for c = 0.07 and 
c — —0.05 fm, predict basically the same total neutron 
skin thickness with a similar quality and they are close 
to the average ARnp = (0.90±0.15)/-H (-0.03±0.02) fm 
[48, 49] of the experimental data. However, the splitting 
of the neutron skin thickness into a part coming from 
the distance t and another part coming from the sur- 
face width Ai?^p is different in both cases, as we have 
discussed in the previous section. Therefore, it becomes 
clear that the experimental neutron skin thickness data, 
by themselves, may be able to constrain the total value 
but not its partition into a bulk and a surface width con- 
tribution. 

It is generally acknowledged that the value of the bulk 
nuclear symmetry energy coefficient J is about 30-32 
MeV, but there is some uncertainty. To assess the de- 
pendence of the extraction of J/Q on the assumed value 
for the J coefficient (which in the above we have taken 
as 31.6 MeV), we have repeated the fit of Eq. (11) to 
the neutron skin data for J = 35 MeV and for J = 28 
MeV. The results are, respectively, J/Q = 0.642 ± 0.046 
and J/Q = 0.701 ± 0.048 if c = 0.07 fm, and J/Q = 
0.764 ± 0.048 and J/Q = 0.829 ± 0.051 if c = -0.05 fm. 
All in all, it seems safe to consider that within the present 
model the J/Q values compatible with the data from an- 
tiprotonic atoms span a window from about J/Q = 0.59 



to J/Q = 0.88, or 0.6 < J/Q < 0.9 in round figures. 

We recall that the DM of nuclei does not incorpo- 
rate shell effects and averages the corresponding quantal 
magnitudes. With the method described we have fitted 
ARf^P that in general contains shell effects, and possi- 
ble correlation and deformation contributions. But, as 
mentioned, the neutron skin data analyzed show a well 
defined linear trend with the relative neutron excess / 
(namely, ARnp = (0.90 ± 0.15)/ + (-0.03 ± 0.02) fm 
[48, 49]). This trend, and the agreement of the DM val- 
ues for ARnp with the self-consistent ETF calculations in 
finite nuclei, which also are free of shell effects and have a 
linear trend with / (Figs. 2 and 3), gives more reliability 
to the predictions obtained with the DM formula from 
the experimental data. 



B. Constraints on the L parameter 

As we have discussed previously, the parameter L has 
a direct relation with the slope of the symmetry energy 
of the nuclear EOS, see Eq. (8). Having determined the 
values of the J/Q ratio compatible with the experimental 
neutron skins measured in antiprotonic atoms, we can use 
the linear correlation between L and J/Q found in mean- 
field effective nuclear interactions, to obtain an insight on 
the values of the parameter L favored by neutron skins. 

We have displayed the linear correlation L — mJ/Q+n 
in the rightmost panel of Fig. 1 (the linear correlation 
coefficient for the shown interactions is r = 0.978). The 
values of the m and n coefficients have some dependence 
on the set of interactions chosen to make the linear re- 
gression. We have checked that this dependence is rather 
weak. Namely, we have tested the correlation of L with 
J/Q by taking into account successively 10, 14, 18, and 
24 interactions and we have found the linear regressions 
to be comprised between L — 139 J/Q — 52 MeV and 
L = 150 J/Q — 57 MeV. Considering these two limiting 
cases and the constraint 0.6 ^ J/Q 0.9 found in the 
previous section, leads to a variation of L between 31 
MeV and 78 MeV. Thus, our estimate for the L coef- 
ficient, which takes into account the surface width cor- 
rection in ARnp obtained in the calculations with mean- 
field interactions, basically lies in the range 30 < L < 80 
MeV. Had we kept the value of J fixed at 31.6 MeV, 
the extracted range for L would be a little narrower: 
35 < L < 70 MeV. 

In a previous work [51] we have investigated the corre- 
lations between the symmetry energy coefficient in finite 
nuclei and in the EOS at subsaturation densities. These 
correlations allow one to derive an aproximate formula 
for the neutron skin thickness with explicit dependence 
on the L coefficient. By comparison of that result for the 
neutron skin thickness with the experimental data set of 
Refs. [48, 49], in Ref. [51] we found a range of values 
L = 55 ± 25 MeV (displayed in Fig. 3 of Ref. [51]) when 
one includes the surface width contribution Ai?f™ in the 
calculations. That prediction is consistent with the val- 
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ues obtained here by the present procedure. A somewhat 
higher range L = 75 ± 25 MeV was obtained [51] when 
one neglects Ai?*™. Althougli a vanishing AR^ value, 
corresponding to 6„ = bp in the nucleon density distribu- 
tions, is not favored by the mean field interactions (see 
Section III), it cannot be discarded without having more 
experimental evidence on the value of &„ as we have noted 
in Section III (see also Refs. [28, 48]). 

In recent years, considerable advances in probing ex- 
perimentally the density dependence of the symmetry 
energy at subsaturation have been achieved in heavy 
ion collisions (HIC) at intermediate energies. It has 
been found that the symmetry energy can be modelized 
around the saturation density with reasonable good ap- 
proximation by [11-20] 



^sym (p) — J 



(13) 



From Eq. (13), one can estimate the parameter L de- 
fined in Eq. (6) from the stiffness 7 of the symmetry 
energy, as L = 37J. The values for 7 extracted in the 
literature from different HIC observables fall in the range 
7 ^ 0.55 — 1.05, which implies L values roughly between 
50 and 100 MeV. In these studies, the value of J in 
Eq. (13) normally has been taken equal to 31.6 or 32 
MeV. The extraction of the equation of state of cold nu- 
clear matter from HIC data is a very complicated task 
and requires model assumptions [9-21, 69-71]; the indi- 
cated estimates for 7 and L may be somewhat modified 
as more measurements and analyses be performed. The 
range 30 < i < 80 MeV of L values determined here from 
neutron skins, assuming the dependence of Eq. (13), fa- 
vors a constraint 0.32 < 7 < 0.84 for the 7 exponent. 
Thus, the result points towards a soft symmetry energy. 

Our estimates for the stiffness 7 can be compared with 
alternative predictions derived in the recent literature. 
For instance, our range 0.32 < 7 < 0.84 overlaps with 
the values 7 0.55-0.77 that Danielewicz [29] obtains 
from the study of binding energies, neutron skins, and 
isospin analog states of selected nuclei. It also contains 
the value 7 ~ 0.55 that is inferred from the analysis 
of neutron-proton emission ratios in HIC carried out by 
Famiano et al. [20], as well as with the value 7 ^ 0.69 
obtained by Shetty et al. [17-19] from isotopic scaling in 
intermediate-energy nuclear reactions. 

The stiffness of the symmetry energy at subsaturation 
densities also has been investigated from isospin difi'usion 
data in HIC, by means of simulations with an isospin- and 
momentum-dependent transport model with in-medium 
nucleon- nucleon cross sections [12-16]. In this case the 
prediction is 0.69 ^ 7 ^ 1.05. It corresponds to a behav- 
ior of Cgymip) that nearly ranges between the familiar 
p^/^ dependence of the purely kinetic symmetry energy 
of a free Fermi gas, in the lower limit, and linearity in 
the density p, in the upper limit. The constraints on 
the stiffness of the symmetry energy derived from isospin 
diffusion, combined with an analysis of the properties of 
Skyrme interactions, are found to lead to a constraint 
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FIG. 6: Comparison of the estimated values of tlie parame- 
ter L from different observables and methods. Some of the 
estimates have been analyzed through Eq. (13) for Csym{p)- 



63 < L < 113 MeV [12-16]. The predictions from isospin 
diffusion are thus a little stiffer, though the lower limits 
of 7 and L obtained using this method are in agreement 
with the upper limits of 7 and L obtained from our study 
of neutron skins with inclusion of the surface width con- 
tribution. 

Another valuable reference comes from the celebrated 
Thomas- Fermi model of Myers and Swi§,tecki [64, 72]. 
This model was fitted very precisely to the binding en- 
ergies of a comprehensive set of 1654 nuclei. It predicts 
an EOS that leads to a coefficient L = 49.9 MeV. Note 
that if we compare Caym{p) calculated from the EOS of 
the Thomas-Fermi model with Eq. (13), an exponent 
7 = 0.51 is obtained. Additional information on the 
density content of the symmetry energy arises from the 
constraints on the symmetry pressure Psym = PoL/'i ex- 
tracted by Klimkiewicz et al. [24] from the properties of 
pygmy dipole resonances in nuclei. These are indicative 
of a value 7 ~ 0.35-0.65 if one assumes Paym = Pol J 
following from Eq. (13) given above. On the other 
hand, Trippa et al. [25] have obtained the constraint 
23.3 < Csj/to(/9 = O.lfm"^) < 24.9 MeV from consider- 
ation of the giant dipole resonance in ^'^^Pb, which im- 
plies a range ~ 0.5-0.65 for the 7 exponent. We depict in 
Fig. 6 the estimated ranges of values for the L parameter 
from the discussed analyses. 

In summary, in spite of the discrepancies in the details, 
the various findings from experimental isospin-sensitive 
signals, including ours, agree all on a rather soft nu- 
clear symmetry energy at subsaturation densities. Re- 
cent studies of pure neutron matter at low densities based 
on universal properties of dilute Fermi gases lead to a 
similar conclusion [73, 74]. One may mention that there 
exists recent circumstantial evidence [75], derived from 
7r^/7r+ ratios in central HIC collisions at SIS/GSI ener- 
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gies, hinting at that the nuclear symmetry energy is soft 
also in the regime of suprasaturation densities (p > 2po)- 
However, further experimental and theoretical confirma- 
tions of this fact need to be awaited [75]. 



V. SUMMARY AND CONCLUSIONS 

The droplet model predicts that the neutron skin thick- 
ness of atomic nuclei is correlated with the ratio J/Q, 
where J is the symmetry energy in bulk matter and Q 
is the surface stiffness coefficient. We have shown that 
the J/ Q ratio displays a linear relationship with the DM 
parameter L in nuclear mean field models that are cal- 
ibrated to experimental ground-state properties such as 
binding energies, charge radii, and single-particle data. 
In this way, the known Correlation between the neutron 
skin thickness in a heavy nucleus and the density deriva- 
tive of the symmetry energy (or of the neutron equation 
of state) evaluated at a subsaturation density, can be in- 
terpreted in the context of the DM. 

According to the droplet model, the neutron skin thick- 
ness is correlated with the overall relative neutron excess 
/ = {N — Z)/A of nuclei. This fact is in agreement 
with the experimental findings using information from 
antiprotonic atoms. The DM expression for the neutron 
skin thickness contains "bulk" and "surface" parts. The 
bulk part corresponds to the contribution proportional 
to the distance t between the nciitron and proton mean 
surface locations. This part is quite dependent on the 
Skyrme force or RMF parameterization used to compute 
it (see Fig. 2). In finite nuclei, this DM bulk contribution 
systematically underestimates the neutron skin thickness 
extracted directly as the difference of the neutron and 
proton rms radii of the nucleus from self-consistent ETF 
calculations with effective forces. This evidence indicates 
that the surface part, due to the 6„ ^ bp contribution, 
is necessary in the DM formula to properly estimate the 
neutron skin thickness in finite nuclei in mean-field mod- 
els using effective nuclear interactions. 

The DM surface contribution AR^ to the neutron 
skin thickness is smaller than the bulk part and shows 
a well defined linear increasing tendency with the over- 
all relative neutron excess /. We have investigated the 
dependence of the slope u*"' with respect to I of this sur- 
face contribution using various Skyrme and RMF forces. 
We have found that the slopes o"""' lie in a region of 
the a^'^-J/Q plane that can be roughly limited by two 
straight lines as a function of the J/ Q value. It implies 
that nuclear properties, which are included in the cali- 
bration of the free parameters of the Skyrme and RMF 
interactions, constrain the possible values of the surface 
width contribution to the neutron skin thickness in the 
DM. If the same nuclear interaction is used, a good 
agreement between both the DM formula and the self- 
consistent ETF calculations of ARnp in finite nuclei is 
found along the whole periodic table when the contribu- 
tion Ai?^™ is included in the DM. 



To analyze possible bounds suggested by experimen- 
tal neutron skin data on the value of the J/Q ratio, we 
have adjusted the DM neutron skin thickness formula to 
the neutron skin sizes measured in antiprotonic atoms 
[48, 49]. We have determined a window 0.6 < J/Q < 0.9 
for the J/Q ratio by using the largest and smallest surface 
contributions AR^ obtained from successful Skyrme 
forces and RMF parameterizations. These two fits repro- 
duce the experimental data with almost the same qual- 
ity. In other words, the experimental data of the neutron 
skin thickness in finite nuclei constrain the total theo- 
retical estimate but not its partition into a bulk and a 
surface contribution. Once the window of J/Q values is 
known, the compatible range of values of the parameter 
L can be estimated from the linear correlation between 
L and J/Q shown in Fig. 1. From our analysis we find 
the constraints 30 < L < 80 MeV. 

If a model symmetry energy Csymip) = J{p/ Po)'^ is as- 
sumed, a prediction for the value of the stiffness 7 of the 
symmetry energy can be obtained, with use of the empir- 
ical values of J and pQ. In this way we find the estimate 
0.32 ^ 7 < 0.84. Thus, our analysis of the experimen- 
tal neutron skins deduced from antiprotonic atoms sug- 
gests a relatively soft symmetry energy, in good accord 
with the recent indications from pygmy [24] and giant [25] 
dipole resonances. Our prediction for the stiffness 7 of 
the symmetry energy also is in reasonable agreement with 
the constraints derived by Danielewicz [29], by Famiano 
et al. [20], and by Shetty et al. [19] using different obscrv- 
ables, as well as with the value 7 = 0.51 of the EOS of the 
Thomas- Fermi model of Myers and Swi^tecki [64, 72] . In 
the upper Hmit, our prediction overlaps with the lower 
limit provided by the analysis of isospin diffusion data in 
intermediate-energy heavy ion collisions [12-16]. 

In summary, different techniques of extracting the pa- 
rameters that describe the density dependence of the 
symmetry energy predict different values. However, tak- 
ing the average of the central values of the predictions 
displayed in Fig. 6, with account of their uncertainties 
when available, there exist narrow windows of the pa- 
rameters L ~ 45-75 MeV and 7 ~ 0.5-0.8 which are, 
actually, compatible with all the different methods men- 
tioned to obtain them. These ranges of values of the 
indicated parameters for describing the leading density 
dependence of the symmetry energy in bulk matter seem 
to be the "optimal" ones according to present experimen- 
tal evidence from nuclear data. 

To conclude, we are aware that the neutron skin thick- 
ness data derived from antiprotonic atoms are to some 
extent model dependent, and have for some nuclei large 
error-bars. Also, our theoretical method represents just 
an average approximation. In spite of these limitations, 
we hope to have shown that from neutron skin data it is 
possible to make a reasonable estimate of the density de- 
pendence of the symmetry energy with uncertainties that 
are not significantly much larger than those currently ob- 
tained from other experimental observables. One expects 
that future data from the planned parity-violating elec- 
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tron scattering experiment for measuring the neutron ra- 
dius in ^"®Pb [2] will contribute to narrow down the con- 
straints derived here from the thickness of neutron skins 
of nuclei. 
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VI. APPENDIX 

To compute the surface stiffness coefficient Q, and also 
the neutron and proton surface widths 6„ and bp that 
appear in the contribution ARf^ of Eq. (10) to the neu- 
tron skin thickness, we obtain the self-consistent neutron 
and proton density profiles in asymmetric semi-infinite 
nuclear matter ( ASINM) . To do that we consider a semi- 
infinite slab with a plane interface separating a mixture 
of protons and neutrons at the left, whose densities de- 
crease smoothly to zero at the right as empty space 
is reached. The axis perpendicular to the interface is 
taken to be the z axis. Thus, the relative neutron ex- 
cess 6 = {pn — Pp)/{Pn + Pp) depends locally on the z 
coordinate. When z goes to minus infinity, the neutron 
and proton densities approach the values of asymmetric 
uniform nuclear matter in equilibrium, corresponding to 
an interior bulk neutron excess Jo- 
in order to obtain the proton and neutron densities in 
ASINM one has to minimize the total energy per unit 
area with respect to arbitrary variations of the densities, 
with the constraint of conservation of the number of pro- 
tons and neutrons. When do is not very large, so that 
occurrence of drip nucleons does not take place (which is 
the situation in all cases considered in the present work) , 
the constrained energy per unit area reads [56, 58, 68] 



E 
5 



/OO 
[e{z) - llnPn{z) - lJ'ppp{z)]dz. 
-OO 



(14) 



In this equation, e{z) represents the energy density func- 
tional of the nuclear effective interaction under investiga- 
tion, and /u„ and /Xp are the neutron and proton chemical 
potentials. The explicit expressions of s{z) in the ex- 
tended Thomas-Fermi (ETF) method for Skyrme forces 



and relativistic mean field interactions can be found in 
e.g. Appendix A of Ref. [58]. The proton and neutron 
densities obey the coupled local Euler-Lagrange equa- 
tions 



6e{z) 

Spn 



IJ'n = 0, 



Sejz) 
Spp 



0. 



(15) 



We solve them fully self-consistently by numerical itera- 
tion [62]. We note in this respect that we have not used 
any parameterized form of the densities such as e.g. Fermi 
shapes. In the relativistic model the variational equa- 
tions (15) are supplemented with additional field equa- 
tions for the meson fields; the calculational details for the 
relativistic problem can be found in Refs. [58, 62, 63]. 

From the calculated density profiles, one obtains the 
mean locations of the surfaces, zoq [q = n,p), as 



i-ooZp'q{z)dz 
I-^Pqiz)dz ' 



(16) 



where the primes indicate a derivative with respect to the 
z coordinate. The proton and neutron surface widths in 
ASINM are obtained as the second moment of p'{z): 



K = 



I-oo(^ - ZOq)^p'g{z)dz 
IZoP'qiz)dz 



(17) 



The distance t = zon — zop between the mean surface 
locations of the neutron and proton density profiles al- 
lows one to extract the surface stiffness coefiicient Q from 
the relation [52] 



r — zon — zop — ~2~ o 



(18) 



which is valid in the limit of small asymmetries. For each 
given nuclear interaction we solve the ASINM problem 

for 5 different values of Sq, between 0.005 and 0.025, and 
we then evaluate Q from the slope of t. 

There exists a second way of computing Q. It is based 
on the fact that in dividing the energy in bulk and sur- 
face parts, as soon as Sq ^ 0, there arc two possibilities 
to define the bulk reference energy [58, 68]. One defi- 
nition is based on the chemical potentials /x„ and p,p of 
each nucleon species; this definition is thermodynami- 
cally consistent and it is the one that we have given in 
Eq. (14). The second definition of the reference energy 
is based on taking the value of the energy per particle 
in bulk asymmetric nuclear matter. Accordingly, upon 
the bulk reference energy chosen, there exist two forms 
of the surface energy in ASINM, which are called Esuif,ij. 
and -Esurf.o [58, 68]. In the small asymmetry limit, it can 
be shown that the difference between these two quantities 
behaves as 



Thus, the slope of E, 



9J2 2 



(19) 



surf.c - -Esurf.p '^i*'^ rCSpCCt tO 

provides another means to extract the value of the surface 



12 



stiffness coefBcicnt Q. We have computed Q from Eq. our calculated values. 
(18) and have used Eq. (19) to confirm the validity of 



[1] I. Angeli, At. Data Nucl. Data Tables 87, 185 (2004). 
[2] R. Michaels, P. A. Souder, and G. M. Urciuoli, 

spokespersons, Jefferson Laboratory Experiment E06002, 

http : / /hallaweb .j lab .org/ parity /prex. 
[3] J. M. Lattimor and M. Prakash, Science 304, 536 (2004); 

Phys. Rep. 442, 109 (2007). 
[4] H. Hciselberg and M. Hjorth-Jensen, Phys. Rep. 328, 237 

(2000). 

[5] A. W. Steiner, M. Prakash, J. M. Lattimer and P. J. EUis, 

Phys. Rep. 411, 325 (2005). 
[6] T. Sil, M. Centelles, X. Virias, and J. Piekarewicz, Phys. 

Roy. C71, 045502 (2005), and references therein. 
[7] P. Jofre, A. Reisenegger, and R. Fernandez, Phys. Rev. 

Lett. 97, 131102 (2006); P.G. Krastev and Bao-An Li, 

Phys. Rev. C 76, 055804 (2007). 
[8] Dao T. Khoa, W. von Oertzen, and A. A. Ogloblin, Nucl. 

Phys. A602, 98 (1996). 
[9] T. Gaitanos, M. Di Toro, S. Typel, V. Baran, C. Fuchs, 

V. Greco, and H. H. Wolter, Nucl. Phys. A732, 24 

(2004) . 

[10] V. Baran, M. Colonna, V. Greco, and M. Di Toro, Phys. 

Rep. 410, 335 (2005). 
[11] M. B. Tsang et al., Phys. Rev. Lett. 92, 062701 (2004). 
[12] A. W. Steiner and B. A. Li, Phys. Rev. C72, 041601(R) 

(2005) . 

[13] L. W. Chen, C. M. Ko, and B. A. Li, Phys. Rev. Lett. 

94, 032701 (2005). 
[14] L. W. Chen, C. M. Ko, and B. A. Li, Phys. Rev. C72, 

064309 (2005). 

[15] L. W. Chen, C. M. Ko, and B. A. Li, Phys. Rev. C76, 
054316 (2007). 

[16] B. A. Li, L. W. Chen, and C. M. Ko, Phys. Rep. 464, 
113 (2008). 

[17] G. A. Souliotis, D. V. Shetty, A. Koksis, E. Bell, M. 
Jandel, M. Veselsky and S. J. Yennello, Phys. Rev. C73, 
024606 (2006). 

[18] D. V. Shetty, S. J. Yennello, G. A. Souliotis, Phys. Rev. 

C75, 034602 (2007). 
[19] D. V. Shetty, S. J. Yennello, G. A. Souliotis, Phys. Rev. 

C76, 024606 (2007). 
[20] M. A. Famiano ct al., Phys. Rev. Lett. 97, 052701 (2006). 
[21] M. B. Tsang ct al., Phys. Rev. Lett. 102, 122701 (2009). 
[22] J. Piekarewicz, Phys. Rev. C66, 034305 (2002); Phys. 

Rev. C69, 041301(R) (2004). 
[23] T. Li et al., Phys. Rev. Lett. 99, 162503 (2007). 
[24] A. Klimkiewicz et al, Phys. Rev. C76, 051603(R) (2007). 
[25] L. Trippa, G. Colo, and E. Vigezzi, Phys. Rev. C77, 

061304(R) (2008). 
[26] H. Liang, N. Van Giai, and J. Meng, Phys. Rev. Lett. 

101, 122502 (2008). 
[27] W. D. Myers and W. J. Swi§,tecki, Nucl. Phys. A336, 

267 (1980). 

[28] W. .L Swiatecki, A. Trzcihska and J. Jastrz§bski, Phys. 

Rev. C71,' 047301 (2005). 
[29] P. Danielewicz, Nucl. Phys. A727, 233 (2003). 
[30] P. Danielewicz and J. Lee, Nucl. Phys. A818, 36 (2009). 
[31] C. J. Pethick and D. G. Ravenhall, Nucl. Phys. A606, 



173 (1996). 

[32] S. Mizutori, J. Dobaczewski, G. A. Lalazissis, W. 
Nazarewicz and P.-G. Reinhard, Phys. Rev. C61, 044326 

(2000) . 

[33] R. J. Furnstahl, Nucl. Phys. A706, 85 (2002). 

[34] J. Treiner and H. Krivine, Ann. of Phys. 170, 406 (1986). 

[35] M. Warda, B. Norlo-Pomorska and K. Pomorski, Nucl. 

Phys. A635, 484 (1998). 
[36] A. E. L. Dieperink, Y. Dewulf, D. Van Neck, M. Waro- 

quier, and V. Rodin, Phys. Rev. C68, 064307 (2003). 
[37] B. A. Brown, Phys. Rev. Lett. 85, 5296 (2000). 
[38] S. Typel and B. A. Brown, Phys. Rev. C64, 027302 

(2001) . 

[39] M. Centelles, M. Del Estal, X. Vifias, and S. K. Pa- 
tra, in The Nuclear Many-Body Problem 2001, Vol. 53 
of NATO Advanced Studies Institute Series B: Physics, 
edited by W. Nazarewicz and D. Vretenar (Kluwer, Dor- 
drecht, 2002), p. 97. 

[40] M. Baldo, C. Maicron, P. Schuck and X. Vinas, Nucl. 
Phys. A736, 241 (2004). 

[41] S. S. Avancini, J. R. Marinelli, D. P. Menezes, M. M. W. 
Moraes, and C. Providencia, Phys. Rev. C75, 055805 
(2007). 

[42] L. Ray, G. W. Hoffmann, G. S. Blanpied , W. R. Coker, 
and R. P. Liljcstrand, Phys. Rev. C18, 1756 (1978); L. 
Ray and G. W. Hoffmann, Phys. Rev. C31, 538 (1985). 

[43] V. E. Starodubsky and N. M. Hintz, Phys. Rev. C49, 
2118 (1994). 

[44] S. Karataglidis, K. Amos, B. A. Brown and P. K. Deb, 

Phys. Rev. C65, 044306 (2002). 
[45] B. C. Clark, L. J. Kerr and S. Hama, Phys. Rev. C67, 

054605 (2003). 

[46] A. Krasznahorkay et al., Phys. Rev. Lett. 82, 3216 
(1999). 

[47] A. Krasznahorkay et al, Nucl. Phys. A731, 224 (2004). 
[48] A. Trzcihska, J. Jastrzcbski, P. Lubiriski, F. J. Hartmarm, 

R. Schmidt, T. von Egidy and B. Klos, Phys. Rev. Lett. 

87, 082501 (2001). 
[49] J. Jastrz§bski, A. Trzcihska, P. Lubihski, B. Klos, F. J. 

Hartmann, T. von Egidy and S. Wycech, Int. J. Mod. 

Phys. E13, 343 (2004). 
[50] B. Klos ct al., Phys. Rev. C76, 014311 (2007). 
[51] M. Centelles, X. Roca-Maza, X. Vifias, and M. Warda, 

Phys. Rev. Lett. 102, 122502 (2009). 
[52] W. D. Myers and W. J. Swi9,tecki, Ann. of Phys. (N.Y.) 

55, 395 (1969); Ann. of Phys. (N.Y.) 84, 186 (1974). 
[53] W. D. Myers, "Droplet Model of Atomic Nuclei" 

(Plenum, New York, 1977). 
[54] M. Brack, C. Guot and H.-B. Hakansson, Phys. Rep. 123, 

275 (1985). 

[55] M. Farine, J. M. Pearson and B. Rouben, Nucl. Phys. 
A304, 317 (1978); M. Farine, J. Cote and J. M. Pear- 
son Nucl. Phys. A338, 86 (1980); Phys. Rev. C24, 303 

(1981). 

[56] K. Kolehmainen, M. Prakash, J. M. Lattimer, and J. 

Treiner, Nucl. Phys. 439, 535 (1985). 
[57] M. Centelles and X. Vinas, Nucl. Phys. A563, 173 



13 



(1993). 

[58] M. Ccntelles, M. Del Estal and X. Vifias, Nucl. Phys. 

A635, 193 (1998). 
[59] J. Piekarewicz and M. Centelles, Phys. Rev. C79, 054311 

(2009). 

[60] M. Lopez-Quelle, S. Marcos, R. Niembro, A. Bouyssy, 

Nguyen Van Giai, Nucl. Phys. A483, 479 (1988). 
[61] I. Bednarek and R. Manka, arXiv:0905.0131. 
[62] M. Centelles, M. Pi, X. Vinas, F. Garcias and M. Bar- 

ranco, Nucl. Phys. A510, 397 (1990). 
[63] M. Centelles, X. Vinas, M. Barranco and P. Schuck, Ann. 

of Phys. (N.Y.) 221, 165 (1993). 
[64] W. D. Myers and W. J. Swi§,tecki, Nucl. Phys. A601, 

141 (1996). 

[65] J. P. Blaizot, Phys. Rep. 64, 171 (1980). 

[66] S. K. Patra, M. Centelles, X. Viiias, and M. Del Estal, 

Phys. Rev. C65, 044304 (2002). 
[67] G. Colo, N. Van Giai, J. Meyer, K. Bennaceur, and P. 



Bonche, Phys. Rev. C70, 024307 (2004). 
[68] W. D. Myers, W. J. Swiatecki, and C. S. Wong, Nucl. 

Phys. A436, 185 (1985). 
[69] S. Kowalski et al., Phys. Rev. C75, 014601 (2007). 
[70] S. K. Samaddar, J. N. De, X. Vinas, and M. Centelles, 

Phys. Rev. C76, 041602(R) (2007); Phys. Rev. C78, 

034607, (2008). 

[71] S. R. Souza, M. B. Tsang, R. Donangelo, W. G. Lynch, 
and A. W. Steiner, Phys. Rev. C78, 014605 (2008). 

[72] W. D. Myers and W. J. Swi9,tecki, Phys. Rev. C57, 3020 
(1998). 

[73] A. Schwenk and C. J. Pethick, Phys. Rev. Lett. 95, 

160401 (2005). 
[74] J. Piekarewicz, Phys. Rev. C76, 064310 (2007). 
[75] Z. Xiao, B. A. Li, L. W. Clicn, G. C. Yong, and M. Zhang, 

Phys. Rev. Lett. 102, 062502 (2009). 



